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Using non-equilibrium Green's functions, we derive a formula for the electron current through 
a lead-molecule-lead nanojunction where the interactions are not restricted to the central region, 
but are spread throughout the system, including the leads and the lead-molecule interfaces. The 
current expression consists of two sets of terms. The first set corresponds to a generalized Meir and 
Wingreen expression where the leads' self-energies are renormalized by the interactions crossing at 
the molecule-lead contacts. The second set corresponds to inelastic scattering events in the leads 
arising from any arbitrary interaction, including electron-electron and electron-phonon coupling, 
treated beyond mean-field approximations. Using different levels of approximation, we are able to 
recover well-known expressions for the current. We also analyse how practical calculations can be 
performed with our formalism by using the new concept of generalized embedding potentials. 
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I. CONTEXT 

Developing a theory for the non-equilibrium electronic 
quantum transport through nanoscalc junctions is a chal- 
lenging task, especially when thinking in terms of ap- 
plications for nanoscale electronics. Electronic transport 
through nanojunctions (single-molecule junctions, for ex- 
ample) exhibits many important new features in compari- 
son with conduction through macroscopic systems. This 
leads to promising new applications in single-molecule 
electronics. In particular, interactions such as Coulomb 
interactions between the electrons and scattering from 
localized atomic vibrations are critically important. 

Having a simple expression for the current (or the 
conductance) of a nanoscale object connected to ter- 
minals is most useful. This is provided by the Lan- 
dauer formulcP which describes the current in terms of 
local properties (transmission coefficients) of a finite cen- 
tral region C and the distribution functions of the elec- 
tron reservoirs connected to this region C. However, 
the original Landauer formulation deals only with non- 
interacting electrons. It has been used with success in 
conjunction with density-functional theory calculations 
for realistic nanoscale system^^l^ since DFT maps the 
many-electron interacting system onto an effective single- 
particle problem. However there are many cases w hen 
such a single-particle approach becomes questionable^. 
The Landauer formula has been built upon by Meir and 
WingreerP to extend the formalism to a central scat- 
tering region containing interactions by using the non- 
equilibrium Green's functions formalism. Other general- 
izations of Landauer-like approaches to include interac- 
tions and inelastic scattering in the region C have been 
developed^"^^'^. However, in real systems the interaction 
is not confined to the central region but exists through- 
out the system. Acc ounting for the interaction along the 
whole system is vita F * ^^ *^. 

In this paper, we provide a complete description which 
generalizes the Meir and Wingreen formalism to systems 



where interactions exist throughout the system, as well 
as at the interfaces between the central region and the 
electrodes. Since the choice of the location of these in- 
terfaces is purely arbitrary, and since the interactions 
exist everywhere, ou r app roach is formally identical to a 
partition-free scheme^^Eni. While keeping the approach 
of the original work of Meir and Wingreert^, we derive 
the most general expression of the current for the fully 
interacting system. From this, we recover all previously 
derived transport expressions or corrections when intro- 
ducing the appropriate level of approximation for the in- 
teraction. Our formalism also leads to the generalization 
of the concept of embedding potentials when the interac- 
tion crosses at the boundaries. It therefore provides an 
alternative way of introducing open boundary conditions 
with interaction in finite-size systems. 

The paper is organized as follows. In Section [Hj we 
provide the generalised current formula for fully inter- 
acting systems. We describe our model in Section |II A| 
and derive the current expression in Section |IIB[ with 
full details of the calculations provided in Appendix [C] 
The connections between more conventional results: the 
current at equilibrium, the current formula of Meir and 
Wingreen and others are given in Sections II C to IIP 



In section II G we describe how to apply our formalism 
in a specific case of interaction crossing at the contacts. 
Finally, we conclude and discuss extension of our work 
in Section Hill 



II. NON-EQUILIBRIUM QUANTUM 
TRANSPORT 

A. The model 

The system consists of two electrodes, left L and right 
R, which connect a central region C via coupling matrix 
elements. The interaction, which we specifically leave 
undefined (e.g. electron-electron or electron-phonon), is 
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FIG. 1; Schematic representation of a central scattering re- 
gion C connected to the left L and right R electrodes, with 
respective quantum-state labels {A},{n},{p} for the three 
L, C, R subspaces. Interactions are given by the coupling of 
the region C to the L{R) electrode Vlc/cl {Vrc/cr)^ and 
by the many-body effects within all regions as well as 

across the LC and CR interfaces. Top: interfaces are arbitrar- 
ily placed at the contact between the scatterer (a molecule) 
and the leads. They cannot be considered as being at equilib- 
rium, being each in a region of strong spatial variation of the 
current and potential drop. Bottom; interfaces are now well 
inside the L, R regions (C region is now the so-called extended 
molecule), and are at local quasi-equilibrium. 



From the continuity equation Vj -|- dtn = 0, we write 
the current through the interface between the L and C 
regions. The current flowing at the LC interface is given 
by I^it) = -edt{NL{t)), where (iVi(t)) is the number 
of electrons in the L region, and is related to the lesser 
Green's function &sjjhjt)) = J2x -iG^;^(t,i). Fr om the 
equations of motior!^-Sli?' obeyed by the Green's functions 
on the Keldysh time-loop contour Ck, we obtain the cur- 
rent II as 

hit) = ^Tr;, [(SG)<(i,t) - (GI])<(i,t)] (1) 

From the rules of analytical continuation on Ck (see Ap- 
pendix p|, we find that 



(2) 



There are no lesser (greater) components for Vlc since 
its time dependence is local VLc{t,t') — yLc{t)S{t — t'). 

In the steady state, all double-time quantities X{t,t') 
depend only on the time difference X{t — t'). The steady 
state current is given after Fourier transform by using 
(SG)<(i,t)^- /da;/27rE(a;)G(a;). To obtain the current, 
we need to calculate the following trace: 



spread over the entire system and crosses at the interfaces 
between the L{R) and G regions. We use different labels 
for the quantum states on each side of these interfaces: 
{A, A'}, {n, m}, {p, p'} are used to represent the complete 
and orthogonal set of states for the L, G and R regions 
respectively. We also use a compact notation for the 
matrix elements M of Green's functions {g,C), the self- 
energies (S) and coupling to the leads (V), where Mc 
represents the matrix elements Mnm in the region G, 
Mlc for Mxm, McL for M„a', Mrc for Mp^, Mcr for 
M„p,, Ml for Maa', and Mr for Mpp,. 

The complete derivation of the current expressions for 
the fully interacting lead-central region-lead junctions re- 
lies on only two assumptions: the many-body effects 
of the interacting particles are well described by self- 
energies in the one-particle Green's functions G, 
and there is no direct coupling or interaction between 
the states of the L and R regions: the only interaction 
between the leads is mediated by the region G, there is 



no direct coupling i.e. S^^^,^,^ 



(P. 



Ti-A [...]= E ^A„G<,(c.) - G<„(^)Ka 

A,n,7 

+ ^■<(c.)G^a(c.) + E^^«-'^(..)G<,(c.) 



(3) 



where 7 runs only on the L and G regions, since = 
(there is no direct coupling between the L and R regions) . 
We then need to evaluate the Green's functions matrix 
elements G<^,C<^,C'^^ and G^^, and Cf^^/" by using 
the Dyson equation Cfj = gfj + [g'EC]fj (with x — r,a,< 
and the indices for the corresponding matrix ele- 

ments) and the rules of analytical continuation for the 
products (see Appendices [b] and [C] for detail). 

We find the following general expression for the current 
II flowing through the left interface: 



Tr 
+TrA 



duj 
2^ 



Ghfic + GUric)^ + G<(tic - tl^)l (4) 



j-,MB,>^< 



A'A 



Y^MB.<^> 



B. The non-equilibrium current 

The location of the interfaces LC and CR is arbi- 
trary (Fig. [1]) , and chosen purely for mathematical con- 
venience, as the interaction spreads throughout the sys- 
tem. We include such interfaces to make the connection 
between our results and other previously derived expres- 
sions within the partitioning scheme. 



where 



'Y't v*^ Pi"" 

^ LC ~ ^CL tJL ^LC^ 

"^LC = ^CL (ffi ~ 9l) ^LC + ^CL 9l ^lc- 



(5) 



By definition 'E^ci'-^) = ^lc + ^^lc (i^): similarly for 
the CL components. g^£°'{uj) are the Green's functions 
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of the region L renormalised by the interaction inside 



where all 



that region: {gT)-^ = H^'')-^ - E 
quantities are defined only in the subspace L. 

There are two contributions to 1^: the first trace is a 
generalisation of the Meir and Wingreen expressiorP to 
the cases where the interactions exist within the three 
L, C, R regions as well as in between the regions. The 
different quantities T^c are related to the generalised 
embedding potentials (i.e. lead self-energies) with inter- 
action crossing at the LC and CR interfaces (see end of 
Appendix [C]) . The second trace in Eq. Q is related to 
inelastic effects involving a sum over the states of the L 
region. Although the L region is semi-infinite by defi- 
nition, an appropriate choice of the location of the LC 
interface reduces the summation. For a closed system at 
equilibrium, the trace Tr [S^G^ — E^G^] is zero simply 
because the system obeys the detailed balance equation: 
E>G< = E<G>. For all other conditions, if the LC in- 
terface is located deep enough in the L electrode, the sys- 
tem is locally at quasi-equilibrium, and hence the trace 
vanishes (see below). 

An expression similar to Eq. Q can be obtained for 
the current Ir flowing at the right CR interface by swap- 
ping the index L R and using the current conservation 
condition II + Ir — 0- 



Different Green's functions 

Finally we need to know, for practical calculations, the 
different Green's functions in all three regions. To eval- 
uate the currents Il,r, we need the Green's functions 

C(j and . 



We find for Gj; = (n|G'^|m) 



ci^) ^ 9c+9c 



9C ^L+R 



L+R 



{CO) 



(6) 



where i^+i? ^^le sum of the generalised leads' self- 
energies [a = L, R) defined as Y^ = {J^caga^acY ■ 



We find for G^j = (pIG^/o') that 



C^ _ -r I C^ 
^ pp' ilpp' ' ilppi C,piP2 P2 



(7) 



where Y^ is the embedding potential arising from the 
central region C. It is defined in the right region R as 
follows: 



(8) 



with = Vrc + ^^c'^ (similarly for EJ^^^). is a 
retarded Green's function of the region G renormalized 
by the interaction inside the central region G and by the 
embedding potential of the left region L only: 



~fc= [rc{^))-'-Yl{u) 



The form of these equations hold for the Green's func- 
tion G£ as well as for the advanced Green's functions 



Ga 
L,R- 



We finally get for G^ 



(n|G<|TO): 



G^ 



CI 



nMB,< 



Y< 

^ L+R 



Ga 



(10) 



with Y^_^j^iuj) = J2a=L,R(^Caga^ac)^- The rulcs of 
analytical continuation need to be applied to the prod- 
ucts {'Ecai^)gai^)^aci'-^))^'^'''^ to get the full expansion 
of the generalised embedding potentials. 



C. The current at equilibrium 

One of the obvious checks to perform is that there is no 
net current at the LC and CR interfaces at equilibrium. 
Considering the equation for 1^ given by Eq.Q, we have 
already shown that the trace TrA[...] vanishes at equilib- 
rium because of the detailed balance principle. Now we 
have to prove the same for the trace Tr„[...] in Eq.Q. For 
this we use the procedure which consists of introducing 
non-equilibrium distribution functions (see Section II F 
below). Since at equilibrium all distributions are equal 
to the Fermi distribution we end up, after long but 
trivial manipulation of Eq. Q , with 



Tr, 



leq 



Tr„ [G^E^^ {n~gl~~9in)^lc 

-c^c^^cL {n~gi-~9in)^ic], 



(11) 



which after further manipulation (using complex- 
conjugate relations between Green's functions and self- 
energies ) can be shown to be equal to zero. Hence, as 
expected, the current from Eq.Q vanishes at equilib- 



D. Recovering the Meir and Wingreen current 
formula 



For systems where there are interactions only within 
G, we have E^f^ ^ q E^/j/^^ = Vlc/cl- Then 

9a = ga, and Tlc = Vcl 9l Vlc, '^lc = ^cl gt Vlc = 
"C^Lc)^' and Eq. ^ can be recast as 

= 1 1 ^ Tr„ [/i(G^ - CI)Vl + C<Tl\ , (12) 



LC- 



with iJl^l = Vcl 5l Vlc and iF^ = VcL{gl 
Hence we recover the result of Meir and WingreerP. 

Going one step further, we consider interaction within 
the L and R regions as well. The current II in Eq.Q 
tak es th en the form of the Meir and Wingreen expression 
Eq.(12), with renormalised escape rates F^, i.e. iF^ = 

Vcl Cgl-gl) Vlc and = ifLTUio) = Vcl ~gt Vlc- 
The interactions within the leads renormalise the cou- 
(9) pling at the contacts F^. Note that we have allowed for 
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a renormalised distribution function in the definition 
of T^p. The distribution of the left lead /l has the same 
form as the Fermi distribution function, but depending 
on the approximation chosen for the interaction , the 
corresponding Fermi level may also need renormalization. 



E. Transport with interaction on the (TD)DFT 
level 



We consider cases where the interaction is spread 
throughout the entire system, and are treated at the 
level of density-functional theory (DFT). The exchange 
and correlation effects for interacting electrons are given 
by an effective potential Vxc{T,t) obtained from an xc 
action functional of the electron density. To this po- 
tential corresp onds an effective self-energy, local in both 
space and tim^3Il!_ This forms a class of self-energies, 
where E^^(r, r') = S^^(r)(5(r - r') cannot have lesser 
or greater components, since the times r and r' must 
be on the same time-loop branch. With no lesser and 
greater components for E'^^, the trace TrA[...] in Eq.Q 
simply vanishes. We are thus left with 



II = 



2tt 



Tr, 



Gr 'Y'l 
C '^LC 



LC) 



(13) 



where T 

'V"t 



LC 



'^CL 5l ^lc , (Tf 
ScL gl ^LC, and S 



I 

LC 



LC) 
V - 



'^CL gL ^LC, 

c {V has only 
^aC/Ca components, and v^c has local static or dynamic 
components w^c,{A,n,p} for DFT or time-dependent DFT 
calculations respectively JPS! Hence we recover a Meir- 
and-Wingreen-like expression for the current with renor- 
malised Fi. The potential Vxc is spread throughout the 
system and inside the leadJ^l^^. Hence Eq. (13) formally 
confirms the necessity of including the potential drop due 
to Vxc in the linear-response regimd^. 

One should note that our formalism includes all other 
cases with other kind of interactio ns (ele ctron-phonon) 
confined only in the central regiorP^JUZl) jt also in- 
cludes other kind of electron-hole excitations which can 
be present in the leadJ^Sl and provides a way to treat sys- 
tems with electron-hole excitations crossing at the con- 
tacts between the central region and the leads. 



F. The current in terms of distribution functions 
and spectral densities 

We now discuss when the second trace in Eq. Q 
vanishes. We introduce the non-equilibrium distribu- 
tions f^(oj) obtained from the generalised Kadanoff- 
BaymansatJ29!x<(Lj) = /<(a;)X°(w)-X'-(a;)/<(w) for 



a Green's function or a self-energy X. We define 

~r fO< 



'cJc ' 



~gt 


^fl^gl-i 


G<c 


= fcGc ^ i 


^LC 


rint< v^a 

— Jl ^lc 


^CL 


/'int< Y^a 
— Jc ^CL 



.fint< 
^LCJC ' 
•sr^r .fint< 



(14) 



(with the Fermi distribution of the L region) . We 
rewrite Eq. Q as follows 



dw 



Tr 



[5G^ ^CL 



[gl 



~gl) 



o_ ^^ 'f- L — -^O-L vi'lv .'^i^j ^LC\ 

TfA [{^IcGh^^L ~ '^lcGc^cl) ' 
(27r)2Tr, [5f< Afiu^)Afiu;)] , 



~gt] 



(15) 



^< - Sf< G'^r 



^Gl 



c ^c ^fc ■ 

fO< _ Ant< 
JL J L ' 



with<5.9<=<5/^ -gl^-gl Sfl^JG^^ 
The differences of distributions are 5/°*^ 
^ftc ~ fLC ~ /l"c^' ^^'^ spectral functions are 

At equilibrium all distributions are equal to the Fermi 
distribution, all Sf = and = as expected. For in- 



teractions localised in C only, we again recover Eq.( 12 1 by 

VLCiCL) , Af^O and gL 



noticing that S^^T^^^) 



gL- 



Furthermore, when the LC interface is located well inside 
the L region, the states A on the left side of the inter- 
face are at their local equilibrium. Hence the correspond- 
ing distributions are equal to the local Fermi distribution 
and (5/ °^ — Sf^ — 0; and thereforce the traces Tr;^[...] 
in Eq. ^ and in Eq. Q vanish (QED). 
The current expression reduces then to 



I^^lJ ^ Tr„ [SG< {gl " K) ^Ic] (16) 
which is just another way to express the Tr„[...] in Eq.Q. 

G. An example of crossing interaction 



We now give a brief description of how to implement 
our formalism for a specific case. We consider a single- 
molecule junction in the presence of electron-vibron in- 
teraction inside the central region and crossing at one of 
the contacts. We use the following Hamiltonian for the 
central region 



Hr 



{a) + a)dU (17) 



where one electronic level Eq and one vibration mode of 
energy wq are coupled together via the coupling constant 
7o. The central region is coupled to the non-interacting 
L and R regions via hopping integrals Uw- 



Vlc + Vcr= Y1 toaicid + d^ca 



(18) 



OL=L,R 



We also consider that the hopping of an electron from 
the C to the L region (and vice versa) can excite another 
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vibration mode of energy uja via the coupling constant 

Hlc = lA{b^ + b){c{d + d)cL) + ujAb^b. (19) 

This model can be understood as a lowest-order expan- 
sion of the hopping integral toL{X) = toL + t'oL^ be- 
tween the C and L regions in terms of the relative po- 
sition X = yjh/ {2mA'jJA){b'^ + b) of the region C with 
respect to the region L. The Hamiltonian H^c repre- 
sents in this model the interaction crossing at the LC 
interface. The corresponding non-equilibrium Green's 
functions and self-energies can be calculated at differ- 
ent orders of the int eract ion using conventional non- 
equilibrium techniquep£l221^ 

Since there is no other interaction inside the L and R 
regions, the current expression is given by the first line of 
Eq.Q. We consider a mean- field approach to treat the 
crossing interaction. This model leads to the Hartree-like 
expressions for the self-energies at the LC interface: 

^^r^ = -2^i / gG<,(.) (20) 

(similarly for E^]^''^^'^ oc / AujG'^^{uj)). 

The closed expression for Green's function matrix 
elements are calculated from the corresponding Dyson 
equations = [gEG]^^. There are no lesser and 

greater components for the self-energy S^,? at the mean- 
field level, as we have explained in Section |II E| Hence 
Eq.([5]) reduces to 

^ LC ~ ^CL 9l ^lc /21^ 

^ LC — ^CL yL ^LC^ 
C^Lc)^ — "^CL 9l ^LCi 

With L^'p — Eql + ^LC 

One can see that the interaction crossing at the LC in- 
terface induces a static (however bias-dependent) renor- 
malisation of the nominal coupling ipL between the L 
and C regions. This non-equilibrium renormalisation in- 
duces bias-dependent modifications of the broadening of 
the spectral features of the G region. It can also lead to 
new physical non-equilibrium effects in the currenlP^I. 

The effects of the crossing interaction can also be 
treated beyond the mean-fiel d leve l by considering a 
Fock-like dynamical self-energjPsHIl, general the ef- 
fects of other interaction (electron-electron) crossing at 
the LC and/or CL interfaces can be treated in a similar 
manner. 



III. DISCUSSION AND CONCLUSION 

We have derived a exact expression for the current 
through systems with interaction both within the L,C,R 



regions and at the LC and CR interfaces. Our result, 
Eq. ^ , is general, assuming that there are no direct inter- 
actions between the leads; a condition that is physically 
sound, especially for single-molecule junctions where the 
spatial gap between the two electrodes is large enough. 
The location of the LC and CR interfaces with respect to 
the physical realistic scatterer is arbitrary but, in prac- 
tice, should be chosen conveniently for numerical calcula- 
tions. When local quasi-equilibria are reached at the in- 
terfaces, a simpler expression for the current is obtained, 
since the local non-equilibrium distribution functions are 
equal to the corresponding Fermi distributions. The de- 
viations <5/°< and 5f^ represent a quantitative tool to 
determine how far inside the leads the LC/RC interfaces 
need to be to reach local equilibrium. Our formalism pro- 
vides a formal justification of the concept of the extended 
molecule that is commonly used with the conventional 
partitioned scheme. It also provides the correction terms 
needed to deal with interaction crossing at the contacts 
and when the contacts are not in their respective local 
(quasi) equilibrium. 

In practice, the calculations should be performed self- 
consistently since the various self-energies I]'^^ in the 
three regions and at the interfaces are functionals of the 
all Green's functions in all the system. This offers extra 
degrees of freedom to perform non-fully self-consistent 
calculations, and test different levels of approximations 
for the interaction. We have given an example of how 
such calculations can be performed for a specific case in 
Section III Gl We have also found that the current con- 
servation conditions lead to an important result for a 
fully-interacting system: a condition that the many-body 
self-energies Y,^^ should satisfying in order to keep the 
conservation II + Ir = (see Appendix |d| . 

In a broader context, our formalism introduces in a 
formal manner the concept of generalised embedding po- 
tentials to interacting cases. Embedding methods pro- 
vide the correct boundary conditions for solving the 
Schrodinger equation in a limited region of space, region 
I, automatically matching the solution on to the wave- 
function in the rest of the system, region II, via the use 
of the embedding potentiaP^. In the quantum transport 
community, the embedding potentials usually arise from 
the left and right leads to which the scatterer of inter- 
est is connected, and are commonly referred to as the 
lead self-energies. In conventional non-interacting ap- 
proaches, they are given by Sq(w) = Vcaga{^)Vac- In 
our formalism, when the interactions cross the LC/CR 
interfaces, we obtain a generalisation of the embedding 
potentials, defined as Y^{u}) = {T,ca{i^)ga{^)^ac{^))^ ■ 
These generalised embedding potentials contain a dou- 
ble non-locality, in the sense that the many-body part 
of Sqc ha-s a spatial extent different from that of the 
coupling matrix elements Vac- Hence Ya defines a buffer 
zone, contained between two surfaces whose separation is 
related to the characteristic spatial range of the interac- 
tion self-energy T,^^ = S^^^dx^ - x„|). The generalised 
embedding potential provides a new alternative for intro- 
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ducing open boundary conditions with interaction within 
many-body finite size systems. 

In summary therefore, we have introduced a new for- 
malism for an accurate expression for the electron cur- 
rent in fully interacting systems. The expression is gen- 
eral and takes into account the fact that the interaction 
is crossing through the interface on which the current- 
density is integrated. Numerical implementations of our 
formalism will enable us to study cases in which the long- 
range Coulomb interaction is not sufficiently screened 
between the central region and the electrodes to be ne- 
glected or approximated, or cases in which vibration ex- 
citations at the contacts play an important role in the 
transport properties. 



Appendix A: Relationship and symmetry on the 
Keldysh contour 

The relations between the different components of the 
Green's functions and self-energies on the Keldysh time- 
loop contour Ck are given by: 





= x+- 


^-x+- 


= x-+ 


- X— 




= x+- 


^ -X-+ 


= x+- 






x+- 


' + x- 


= x+- 


+ X-+ 




x^- 


^ -x+- 


= X'' - 





(Al) 



with X''i''2(12) EE G''i''2(12) or I]''i''^(12), and where 
(i = 1,2) is the composite index for space-time loca- 
tion {xi,ti) and rji is the index of the Keldysh time- loop 
contour Ck branch (-1- forward time arrow, — backward 
time arrow) on which the time U is located. The usual 
lesser and greater projections are defined respectively as 
X^ = X^ and X-^ = X ^ , and the usual time-ordered 
(anti-time-ordered) as X* = X++ {X* ^ X — ). 

By definition, complex conjugation of the different 
Green's functions follows the rules: 

G"(l,2) = (G''(2,l))* 
G^(1,2) = -(g^(2,1)^ 

Similar expressions hold for the self-energies 11^(1,2) . 

Appendix B: Rules for analytical continuation 

The rules for analytical continuation from Ck to nor- 
mal real-time make that the following products P(j)(r, r') 
on the time-loop contour, 



P(2) = I AB 

Ck 



P(3) - / ABC 

JCk 

P{n) = / AiA2...An, 
JCk 



(Bl) 



have the following components P^-^ {t, t') on the real-time 
axis (a; = r, a, >, <) 



^1) = J^A^'B^ +A^B'' 



= j^A^B^C + A'B^C + A^'B'^C^ (B2) 

pr _ Ar Ar Ar pa _ Aa Aa Aa 

^ J 1 ^---^n ^(n) ~ J ^1^2---^n- 



Appendix C: Derivation of the current II 

In this appendix, we provide the details of the deriva- 
tion of the main results of this paper, mainly Eq. Q and 
Eq. ([5|. In the following, we choose to use the symbol 
J] for summations in order to have a better graphical 
distinction between the sum-signs and the self-energies 
S. 

From Eq. ([T]) and Eq. we need to calculate the 
following traces: 



TiA [(SG)<] = 

5:((V^.„-HSr)G„A)^+E(^^A'^G 



A' A) 



(CI) 



A, A' 



and similarly 
TrA [(GE)<] 



(GA„(s:ff + Ka))< +E(GavS^I?)< (C2) 



A,n A, A' 

since = and Y.^^ = 0. 

We first consider the sums y 

Yi^fyGyx)^ -^(GaA'S 



MBn< 
X'X) 



A. A' 



A, A' 



^AA' "-^A'A + ^AA' '-'A'A ~ '-'AA'^A'A ^'-'AA'^A'A 



A, A' 



'A'A ^ ^A'A^ 



A, A' 

=5:sf,^'<(G<-G>)vA + (S 

A, A' 

Ev^MB,>^< y^MB,<^> 
^AA' *-^A'A ^ ^AA' *-^A' 

A, A' 

=TrA [yT'^G< - Sr'^G> 



MB/ 
AA' 



MB,> ipMB.<\ /-^< 

-2. jAA'G;^,;^ 



(G3) 



In the first line of Eq. (G3), we used the rules of an- 



alytical continuation. In the second, we have used the 
equivalent of cyclic permutation in the calculation of a 
trace, i.e. swapping the index A and A' in the last two 
terms. This is possible here since the sums and all matrix 



7 



elements are defined in the single subspace of the L elec- 
trode. The final result looks like the collision terms usu- 
ally obtained in the derivation of a generalised Boltzmann 
equation from quantum kinetic theory. They correspond 
to the particle production (scattering-in) and absorption 
or hole production (scattering-out) related to inelastic 
processes (i.e. non-diagonal elements of the self-energy 
on the time- loop contour S^) occuring in the left elec- 
trode. 

Now we consider the sums £a n Eq.(Cl| and 
Eq.(IC2l). We find that 



explained above, we have used the condition S 



pA 



0. 



The same principle holds for the derivation of the other 
Green's functions matrix elements. 

The interesting point is that the terms in can be 

factorised out and included within the renormalisation of 
the left lead Green's functions gx{\^ as follows 



9\x' 



(1 - EMB,a^aN 



A'Ai - 5aAi- 



(C6) 



An 
An 



'^An'-^nA 



G< 
Xn^nX 



^MB< 
-'Xn 



nX ^ <-f 



MB<1 



An^riA 



(C4) 



We now need to calculate the following different 
Green's functions matrix elements G5, , Gfv, , Gt >, and 



G\^. For this we use the Dyson- like equation defined 
for the non-diagonal elements: G^y^ = {n\{G'Eg)^\X) , and 
= (A|(GS5)1"), and 
We concentrate on one matrix element {n\{GI]g)'^\X) 
to show the mechanism of the derivation: 



G<^^{n\iGJ:g)<\\) 

= (n|G''S<.g'' + G<SV + G'^E''g<|A) 

^nXi ^AiA2 5a2A + ^mri ^mAg ^AaA 

Ai , A2 ,m 

+ '-^nAi ^AiAa ffAaA + ^nm ^mAa ^AaA 
+ GJ^Ai ^AiAa ffAaA + ^nm ^niAa ^AaA' 



(C5) 



with E 



a/r 



^mA ^^'^ ^mX 



^mA'^' and, as 



Therefore the matrix G^^^ = (n|(GEg)<|A) can be re- 
cast as G'^y^ = {n\(Gc ^cl 5l)^|A), or similarly with an 
explicit summation: 



^nX ^nrn ^mA' 9x'> 



7n.X' 



(C7) 



+ '^nm ^mX' 9X'X + "-^nm ^mAa 9X'X- 

We also find that 



^Xn 

Ga 
nX 



(A|(5L Sic Gcr\n), 
{n\{Gc ^CL ffL)1A), 
GL = (A|(5L GcYW). 



(C8) 



Using the rules of analytical continuation for products 



of three quantities, we find that Eq.(C4) becomes 



Y^{\\T.lc\n) {n\{Gc T^cl ffL)<|A) - (A|(5l ^lc Gc)<\n) (n|E^JA) + (A|Ef<5<|«) (n|(Gc ^cl 5l)1A) 
-(A|(5L Slc GcY\n) (n|I]^f<|A) 

(G< .92 + G^ gl + ~g<) ^IcW) + (n|G^ ^tc " ^Ic Gh\n) 

n 

+ {n\j:^cL {~9l ^Ic + ~9l ^tc + 9l ^Ic G^) \n) 

n 

+ G< (E^i 52 - S2c ~gl S2c) l«) 



Tr„ 



(C9) 
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In the second equality of Eq. ( C9 ) , the matrix elements 



have been swapped to get a trace only over the central 
region subspace {n}. The final two equalities are just 
exactly the Ti\, entering the definition of the current 
given by Eq.(U| with the definitions of T^^^ and Tlc 
(and their adjoints) given by 



^ LC — ^CL 



JCL ^LC 



m 



and 



(CIO) 



(Cll) 



The first term represents the initial conditions . 

For the central region, we have chosen the initial condi- 
tion such as {n\g^\m) = (see Eq.(lOl). We could have 
chosen another initial condition. Such choices have no 
effects on the steady state regime when a steady current 
fiow through the central region, however the initial con- 
ditions play an important role in the transient behaviour 
of the curreni!22H2l!_ 

For the definition of the lesser left and right Green's 
functions, it is however not possible to neglect the initial 
conditions (before full interactions and coupling to the 
region central are taken into account). This is because it 
would not be physical to ignore the presence of the left 
and right Fermi seas, obtained as the thermodynamical 
limit of the two semi-infinite leads which act as electron 
emitter and collector in our model device. 



One can however recast Eq.(Dl ) as follows 



T 



LC 



(QED) Eqs. ( |C3|C9|C10|C11[ ) are just the main results 
of this paper. 

Now we can follow the same mechanism of derivation 
to obtain an expression similar to Eq. Q for the current 
Ijj flowing at the right CR interface. In concrete the 
expression for is given by Eq. Q by swapping the 
index L ^ R and with a minus sign because of the current 
conservation condition II + Ir — 0. 

Finally one should note that because of the following 
three conditions: (i) the very existence of the interaction 
crossing at the contact, (ii) the fact that in the most gen- 
eral cases ^Ca/ac 7^ ^Ca/aC opposition with the non- 



interaction case where / ^ 



V7 



Ca/aC — ^Ca/aC, 

(Hi) the rules of analytical continuation for triple prod- 
ucts P(3) , the usual cyclic permutation performed in the 
calculation of the trace Trx[{TiG)'^ — (GT,)^] cannot be 
used to transform the initial trace over {A} onto a trace 
over {n}. Therefore the current at the LC contact is 
not given by a straightforward generalisation of the Meir 
and Wingreen formula of the type 



and 



duj 
2^ 



Tr„[r<G> 



Tr,[Sr>G< 



(C12) 



where Y£ is the generalised (interacting) embedding po- 
tential of the L electrode. 

This is another very important result of our work which 
has strong implication in the expression of the current 
itself, and also in the conditions of current conservation. 



Appendix D: Current conservation condition 

First we consider the general definition of the lesser 
and greater Green's functions : 

G^ = {1 + G''S'')g^(l + S° G") -I- G'' G^ (Dl) 



G< = G'' ((/)-ig<(5°)-i + I]<) G'^ = G" S< G" (D2) 

with S< = I]< +7< and 7< = (.g'')"\9^(.9")"^- And 
similarly for G^. Hence 'y< ^ j> — (5°)^^ — {9^)~^ and 
t< -t> = (G'^)-! - (G'')-i. 

From these properties, it can be easily shown that 



Trail [S<G>-E>G<] -0 



(D3) 



for each w. The trace runs over all indexes in the system 
(all = {X,n,p}) and the interaction S are spread over 
the whole L, G, R regions. This is the starting point to 
find the conditions for current conservation. 

Because the trace runs over all the three subspaces, 
we can apply the usual cyclic permutation and recast 
Eq.llDSl) as follows 



- Trail [(SG)< - (GE)<] + Traii [^<G> - 7>G<] = 

(D4) 

or equivalently 

J dc^ Trail [(SG)< - (GE)<]+Traii [7>G< -7<G>] = 

(D5) 

Expanding the trace in the first term over each sub- 
space Traii[...] = TrA[..] -t-Tr„[...] -|-Trp[...], one can iden- 
tify the definition of the currents 1^ and //j from Tr;^[...] 
and Trp[...] respectively. 

Hence the condition of current conservation Il+Ir = 
leads to 

J dc.Tr„ [(I]G)< - (GS)<]+Traii [7>G< -7<G>] = 0. 

. 

After further manipulation, lengthy but trivial in the 
light of Appendix [C] we find that the current conserva- 
tion leads to the following condition: 



dLu Tr„^i^a/?,Pr>G<-I]r<G>]-fGic + Ccfl = 

(D7) 
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with 



Clc{^) = Tr„[(f> + Tlc - t[c - T^c)G, 



a.n<iCcR = CLc[{L^R}]. 



(T 



LC 



C 

(D8) 



The first trace in Eq. (D7) corresponds to the sum 



over the three regions L,C,R of the integrated colhsion 
term. The two other traces C^c and Cqr arise from the 
interactions crossing at the LC and CR interfaces. But, 
globahy, Eq. (D7| still implies that the total integrated 



collision terms must vanish. 
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